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ABSTRACT 

We consider gravitational clustering from primoridal non-Gaussian fluctuations provided by 
a x 2 model, as motivated by some models of inflation. The emphasis is in signatures that can be 
used to constrain this type of models from large-scale structure galaxy surveys. Non-Gaussian 
initial conditions provide additional non-linear couplings otherwise forbidden by symmetry that 
cause non-linear gravitational corrections to become important at larger scales than in the Gaus- 
sian case. In fact, the lack of hierarchical scaling in the initial conditions is partially restored 
by gravitational evolution at scales k > 0.1 h/Mpc. However, the bispectrum shows much larger 
amplitude and residual scale dependence not present in evolution from Gaussian initial conditions 
that can be used to test this model against observations. We include the effects of biasing and 
redshift distortions essential to compare this model with galaxy redshift surveys. We also discuss 
the effects of primordial non-Gaussianity on the redshift-space power spectrum and show that it 
changes the shape of the quadrupole to monopole ratio through non-linear corrections to infall 
velocities. 

Subject headings: large-scale structure of universe 
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1. Introduction 

The current paradigm for the formation of large-scale structures in the universe is that small primordial 
fluctuations, with a roughly scale-invariant spectrum, are amplified by non-linear gravitational instability to 
form the galaxies, clusters and superclsuters that we see today in galaxy surveys. A common assumption is 
that the statistics of these primordial fluctuations is Gaussian, which is consistent with current observations. 
However, at this point we cannot rule out the possibility that primordial fluctuations were in fact non- 
Gaussian to some extent. Fortunately, rapid progress in microwave background anisotropy experiments and 
large galaxy redshift surveys in the next few years will provide high quality data that can be used to test 
the nature of primordial fluctuations. 

The consideration of non-Gaussian initial conditions is a difficult issue because there is an infinite class 
of non-Gaussian models. In a Gaussian field <j>, n— point connected correlation functions, £ n = ( (j>\ . . . <f> n ) c of 
order n > 3 vanish; however, in a non-Gaussian field higher-order correlations can virtually behave arbitrarily 
subject to the constraint of realizability, that is, that the multivariate probability distributions be positive 
definite. In general, however, we can divide non-Gaussian models into weakly and strongly non-Gaussian, 
depending on the magnitude of normalized cumulants s p = ( S p ) c / ( S 2 ) p ^ 2 of the smoothed density contrast 
5 compared to unity. A minimal way of achieving strongly non-Gaussian models is by dimensional scaling, 
where the hierarchy of n-point correlation functions obeys oc Q ^ ■ in this case, normalized cumulants are 
numbers not necessarily smaller than unity. This is opossed to weakly non-Gaussian models where normalized 
cumulants are forced to be small quantities, e.g. hierarchical scaling models, where oc Q _1 with ^ 2 <1. 

In this paper we consider \ 2 initial conditions, which belong to the class of dimensional scaling models. 
As a strongly non-Gaussian model, it has the potential to be more easily constrained or detected by large- 
scale structure observations than primordial weakly non-Gaussian models. Furthermore, there is a number of 
inflationary models in the literature that motivate % 2 initial conditions (e.g. Kofman et al. 1989; Antoniadis 
et al. 1997, Linde & Muhanov 1997; Peebles 1997). In addition, this model has been recently argued 
to fit a significant set of observational constraints (Peebles 1999a, 1999b). It is also possible that this 
particular model may be a good representation of the general behavior of dimensional scaling models, and 
thus provide general insight about their advantages and disadvantages. Furthermore, models in which 
primordial fluctuations are generated by topological defects also generally obey dimensional scaling (e.g. 
Turok & Spergel 1991). 

Previous work on clustering from non-Gaussian initial conditions was done using numerical simulations 
(e.g. Moscardini et al. 1991; Weinberg and Cole 1992; Coles et al. 1993; White 1999), and perturbation 
theory (Fry & Scherrcr 1994; Jaffe 1994; Chodorowski & Bouchct 1996; Gaztahaga & Fosalba 1998). In 
this work, we concentrate on aspects of x 2 initial conditions that can be tested with current and future 
observations in galaxy redshift surveys. In particular, we concentrate on the study of the rcdshift-space 
power spectrum and bispectrum; applications of these results to observations will be considered elsewhere 
(Scoccimarro et al. 2000). 

Different aspects of the effects of primordial non-Gaussianity on the power spectrum have been con- 
sidered in the literature (e.g. Feldman, Kaiser & Peacock 1994; Stirling & Peacock 1996; Sutherland et al. 
1999). The skewness of the smoothed density field in texture models was studied by Gaztanaga & Mahoncn 
(1996), whereas the impact of hierarchical scaling models of primordial non-Gaussianity on the bispectrum 
was considered by Verde et al. (1999). Our approach is complementary to recent studies of the impact of 
primordial non-Gaussian models in other aspects of large-scale structure (e.g. Robinson, Gawiser & Silk 
1999; Koyama, Soda & Taruya 1999; Willick 1999; Picrpaoli, Garcia-Bellido & Borgani 1999). 
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This paper is organized as follows. In Section 2 we review perturbation theory results regarding the 
evolution of the bispectrum from general initial conditions. Section 3 presents results for the particular case 
of x 2 initial conditions. In Section 4 we discuss the effects of galaxy biasing and Section 5 those of redshift 
distortions. Finally, we present our conclusions in Section 6. 



2. Gravitational Clustering from Non-Gaussian Initial Conditions 

We are interested in the effects of non-Gaussianity on clustering statistics at large scales. A convenient 
approach is to use non-linear perturbation theory (PT), where the density field in Fourier space at time t 
reads (Fry 1984) 



6{k,t) = 5^{k,t) + 5&{k,t) + ... = D 1 5 I {k) + D\ J S D {k- fe 12 ) F 2 (k 1 ,k 2 )5 I {k 1 )Si{k 2 ) + ..., (1) 

where Si(k) is the initial density contrast, fc^. j = fc, + . . . + kj, and we assumed that to a very good ap- 
proximation, in second-order PT fluctuations grow according to D 2 (t) oc D 2 (t), with Di(t) the linear growth 
factor. The kernel F 2 (k%, k 2 ) = 5/7+1/2 cos8(ki/k 2 +k 2 /ki)+2/7 cos 2 8 with cos 8 = ki-k 2 /(kik 2 ), describes 
to leading order the non-local evolution of the density field due to the long-range nature of gravitational 
interactions. 

The non-Gaussianity of initial conditions is encoded in the statistical properties of the random field 
Si(k), in particular, its lowest order connected moments are 



(*/(*)) = 0, (2) 

(Si(ki)Si(k 2 )) = 6 D (k 12 ) P J (fei), (3) 

(5 I (k 1 )5 I (k 2 )d I (k 3 )) = 6 D (k 123 ) B^kuk^ks), (4) 

(S I (k 1 )S I (k 2 )8 I (k 3 )S I (k i )) c = S D (k 123i ) T / (fc 1 ,fe 2 ,fe 3 ,fc 4 ), (5) 



where P 1 , B 1 and T 1 denote the power spectrum, bispectrum, and trispectrum of the initial density field. 
In linear PT, these just scale with the linear growth factor, but at the scales of current galaxy surveys, non- 
linear corrections can be significant. From Eq. (0) we derive non-linear corrections to the power spectrum 
and bispectrum 

P{k) = P I (k) + 2jd 3 qF 2 (k + q,-q)B I (k,q), (6) 
#123 = B' 123 + Bf 23 + Bl 23 , (7) 

where B\ 23 denotes the contribution of the initial bispectrum scaled to the present time using linear PT, 
^123 W [-Oi(^)] 3 j ^123 represents the usual gravitationally induced bispectrum from Gaussian initial 
conditions (Fry 1984) 



Bg 3 = 2F 2 (k 1 ,k 2 ) P I (k 1 )P I (k 2 )+ eye, 
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and the last term in Eq. (f7j) 

B123 = J A F 2 (k 12 - q, q) T 1 (k l7 k 2 , k 12 -q,q) + eye, (9) 

represents the contribution coming from the initial trispectrum linearly evolved to the present. Note that 
only the first term in Eq. (Q) scales as [-Di(t)] 3 , the last two terms have the same scaling with time, [Di(i)] 4 , 
and therefore eventually dominate at late times. 

3. Evolution from x 2 Initial Conditions 

We now evaluate the results of the previous section for the particular case of % 2 initial conditions. 
In this case, the density field after inflation is proportional to the square of a Gaussian scalar field </>(x), 
p(x) oc 4>(x) 2 . Thus, the density constrast 5(x) = 4>(x) 2 /<t^ — 1, where a 2 , = ( <fi 2 ). If the auxiliary Gaussian 
field has two-point correlation function = ( </>i</>2 ), then the density 2-, 3-, and 4-point correlation functions 
are respectively (Peebles 1999b) 



& = 2^, (10) 



r ?1234 = 4 V ^12^23^34^41 + V ^12^24^43^31 + V Cl3?32^24C41 ; (12) 



/ £l2 £23 £34 £ll + ^^12^24^43^31 + y £l3 £32 £24 £ll 

and the initial density power spectrum, bispectrum, and trispectrum are given by 



P 1 = 2 J d 3 q P^q)P,p(\k-q\), (13) 

B{ 23 = 8jd 3 qP^q)P <t> (\k 1 -q\)P 4> (\k 2 + q\), (14) 
2l234 = 16 J d 3 q Pf(q)Pt(\ki - q\)P4\km - q\)P<t>(\k 3 + q\) + 16 / d 3 q P^q) P^Qk, + q\) 

xP4\k 23 - q\)P^(\k 2 - q\) + 16 J d 3 q P (g)P (|fca + q\)P^(\k 2A - q\)P^(\k 2 - q\), (15) 

where a^P^(k) denotes the power spectrum of the <f> field. For scale-free spectra, P$(k) = Ah™*, 



similarly the bispectrum can be expressed in terms of hypergeometric functions, using the results in Scocci- 
marro (1997). Simple analytic results can be obtained for the particular case — —2, which gives a density 
spectral index n = — 1, reasonably close to the observed one at the non-linear scale. Using Eq. (||) we get 
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However, the initial trispectrum does not seem to have a simple closed form. Defining the non-linear scale 
from the linear power spectrum, A 1 (k n i) = Airk^P 1 \k n i) = 1, we have 



k\ 2 f 24 k\ ( k \ 2 / ._ k 



AW = fcJ l 1 + ^^rteJ l 1 + a77 ^J' <18) 

Rather than working with the bispectrum itself, it is convenient to consider the reduced bispectrum 
Q123 defined by 

Ql23 = 55 wTmTm' (19) 

which for Gaussian initial conditions has the important properties that it is independent of time, and to a 
very good approximation independent of the matter density parameter tt. In addition, for scale- free initial 
conditions is independent of overall scale, for CDM-type models the scale-dependence is only weak through 
the scale variation of the spectral index. From these results, the reduced bispectrum for x 2 initial conditions, 
including non-linear gravitational corrections is 



r , _ 4\/2 k n i 192 k x k 2 + k 2 k 3 + k 3 k x G T 

<^123 - 7 ] 7 —, 7T^ — 71 '. 7 . 7 x 2 rW 123 + yi23) \ Z{J 

tt fei + k 2 + k 3 7ir z (fei + k 2 + k 3 y 

where Q^ 2 3 denotes the reduced bispectrum obtained from Gaussian initial conditions, and Q\ 23 denotes the 
contribution from Eq. (j^) which is difficult to calculate analytically (however, we shall give a full numerical 
evaluation of the bispectrum below). In particular, for equilateral configurations 

g - = 1ST T - S + ? + Q - ra °' 6 X ~ °' 35 + Q - (21) 

where we used that Qf g = 4/7 (Fry 1984). Since the last term is a number independent of scale, Eq. ( pl| ) 
illustrates the signature of this type of non-Gaussian initial conditions: Q eq shows a strong scale dependence 
at large scales as kjk n \ — > 0. This is not just a property of \ 2 initial conditions, but rather of dimensional 
scaling models (£„ oc Q^ 2 )- 

A simple generalization of \ 2 initial conditions is to consider N independent fields, each of them \ 2 
distributed with the same power spectrum P^(fc). As N increases, the Gaussian limit is achieved as a result 
of the central limit theorem. This might be a useful way of parametrizing constraints on non-Gaussianity 
from large-scale structure observations. For a fixed linear density power spectrum P I (k), the primordial 
bispectrum and trispectrum obey B 1 cx 1/y/N and T 1 oc 1/iV, respectively. Equation (^p|) then reduces to 



1 4v/2 k nl 1 192 kxk 2 + k 2 k 3 + k 3 h G 1 T 

!l23 — —7= — — —IT^ — T, — , i \ o r Wi 23 + Wi23- \ zz ) 



N tt h+k-2 + ks N 7~k 2 (fci+fc 2 + fc 3 ) 2 ^ LA6 N 
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Thus, the approach rate to the Gaussian initial conditions result is 1/yN at large N . The same scaling holds 
for the skewness parameter, similarly, the kurtosis relaxes towards the Gaussian initial conditions value as 
S4 oc l/N. Note, however, that the scaling at small N is stronger; in fact, such a behavior is seen in the 
N-body simulation results by White (1999) on the skewness and kurtosis as a function of N (see his Fig. 7). 

From Eq. (|l8|) we see that non-linear corrections to linear PT can be significant even at wavenumbers 
smaller than the non-linear scale. In order to check for significant contributions from third and higher-order 
in the perturbation expansion [Eq. ([!])], we resort to numerical realizations of second-order Lagrangian PT 
(2LPT), which by being formulated in Lagrangian space incorporate the remaining terms in Eq.(Q), although 
only approximately beyond F2 (but the error is small, see e.g. Scoccimarro (1998) for a quantitative com- 
parison of one-point cumulants). In this case, the perturbation expansion is done about particle trayectories 
x(t), 

x(t) = q + ^(q,t) = q + ^ 1 \q,t) + ^ 2 \q,t) + ..., (23) 

so that q represents the initial (Lagrangian) position of a particle whose current (Eulerian) position is x(t) 
and *&(q,t) denotes the displacement vector assumed to be a small quantity. In 2LPT, Eq. ( |23| ) is truncated 
at second order. The solutions for the displacement field are obtained from the equations of motion and read 



x{q) = q-D l V,0W+U a V q {2 \ (24) 
V^ (1) (<?) - S(q), (25) 

V^ (2) (9) = E^^^S^-^S^)) 2 ], (26) 



where to a very good approximation, D2 ~ —3D 2 /7 (Bouchet et al. 1995). For a detailed exposition of 
2LPT see e.g. Buchert et al. (1994) and Bouchet et al. (1995), who also compared to N-body simulations. 
Similarly, in a companion paper (Scoccimarro 2000), we explore the validity of 2LPT for the evolution of 
the bispectrum in redshift-space from Gaussian initial conditions, by comparing to N-body simulations. In 
this paper, we use 2LPT for x 2 initial conditions and adopt the criterion of validity found for the Gaussian 
case, namely, we only include waves up to a maximum wavenumber k c ps 0.5 h/Mpc, so that the percentage 
of shell-crossing is below 10%. Given that 2LPT is computationally inexpensive, we can generate many 
realizations which is essential to beat down cosmic variance. 

Figures |l| and || show the results of 100 2LPT realizations. We have chosen the auxiliary Gaussian field 
</> with a spectral index = —2.4, leading to n = —1.8 as proposed in Peebles (1999a). The amplitude 
of the power spectrum has been chosen to give k n i = 0.33 h/Mpc. First, we checked that the initial 
conditions were generated correctly by comparing the power spectrum and the bispectrum with theoretical 
expectations. The dashed lines in Fig. [I] show the predictions of the first term in Eq. (^0|) for the reduced 
bispectrum at ki = 0.068 h/Mpc, k 2 = 2ki, as a function of angle 8 between ki and fc 2 . This corresponds to 
n = —1, however, it approximately matches the numerical results (triangles, n — —1.8). The latter show less 
dependence on angle, as expected because the scale dependence in the n = —1.8 case (Q 1 oc fc~ - 6 ) is weaker 
than for n = — 1 (Q 1 oc fc -1 ). In Fig. ||we show equilateral configurations as a function of scale for \ 2 initial 
conditions (triangles) and Qi q {k) = 0.8(k/k n i)~ 6 (dashed lines), where the proportionality constant was 
chosen to fit the numerical result, this is slightly larger than the prediction in the first term of Eq. (|l]) for 
n = —1, and closer to the real-space result Q e q{r) = 0.94(r/r n ;) - 6 . 
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The behavior of the x 2 bispectrum is notoriously different from that generated by gravity from Gaussian 
initial conditions for identical power spectrum (dot-dashed lines in Figs. [l|-|]) (Frieman & Gaztahaga 1999). 
The structures generated by squaring a Gaussian field roughly correspond to the underlying Gaussian high- 
peaks which are mostly spherical, thus the reduced bispectrum is approximately flat. In fact, the increase of 
Q 1 as 9 — > 7r seen in Figj^ is basically due to the scale dependence of Q 1 , i.e. as 9 — > tt, the side k% decreases 
and thus Q 1 increases. 

As shown in Eqs. ( po| - |2l| ), non-linear corrections to the bispectrum are significant at the scales of interest, 
so linear extrapolation of the initial bispectrum is insufficient to make comparison with current observations. 
The square symbols in Figs, [l] and ^ show the reduced bispectrum after non- linear corrections are included. As 
a result, the familiar dependence of Q123 on the triangle shape due to the dynamics of large-scale structures is 
recovered (Fig. [I]), and the scale dependence shown by Q 1 is now reduced (Fig. |l]). However, the differences 
between the Gaussian and % 2 case are very obvious: the x 2 evolved bispectrum has an amplitude about 
2-4 times larger than that of an initially Gaussian field with the same power spectrum. Furthermore, the 
X 2 case shows residual scale dependence that reflects the dimensional scaling of the initial conditions. The 
analogous results for the skewness were obtained using numerical simulations (White 1999) and non-linear 
PT in the spherical collapse approximation (Gaztahaga & Fosalba 1998). These signatures can be used to 
test this model against obervations; however, before we can do so we have to test the robustness of these 
conclusions against the effects of galaxy biasing and redshift distortions. 



4. Galaxy Biasing 

We now consider the effects of local biasing when initial conditions are non-Gaussian. If we restrict 
ourselves to scales larger than those relevant to galaxy formation, the galaxy density field can be thought 
of as a local transformation of the density field smoothed over large enough scales so that <5 <§C 1 , and thus 
expanded as (Fry & Gaztahaga 1993) 

S g = b 1 S+ b j-S 2 + ..., (27) 
which implies the following mapping for the power spectrum and bispectrum 

P g (k) = blPW + hbzfdSqB'fcq), (28) 
B g (k u k 2 , k 3 ) = b\B 12i + 6?6 2 S{ 23 + \b\b 2 J d 3 q T( 23 (q), (29) 

where T^io) denotes the trispectrum T I (ki,k2,q) symmetrized over {k±, k 2 , k$}, and X{ 2 3 i s defined in 
Eq. (|l9|). The reduced bispectrum then reads 



b^ 1 " ' b\ 

where the effective non-linear bias parameter is given by 



1 6 off 

Q a = -Qi2Z + ^, (30) 



- 8 - 



LCff 

°2 



or' 
ZZj 123 



d 3 Q Tl 23 (q) 



Qi 



2:S 



^123 



cZ 3 g B 7 (fc 2 ,£?) + cyc. 



(31) 



The analogous result to Eq. ([30|) for the case of the skewness was derived by Fry & Scherrer (1994). So far the 
derivation does not assume any particular non-Gaussian model. In order to compute the effective non-linear 
bias parameter, we need to make assumptions about the initial bispectrum and trispcctrum. The particular 
form of the convolution integrals in Eq. (|^) require knowledge of the three- and four-point functions where 
two points coincide. Using Eqs. (11-12) we obtain 



J d 3 q B'ikx, q) = 2 3 / 2 y/?(p) P'ih), 
J d 3 q T* w (q) = 2 3 / 2 y / e 7 (0) B[ 23 + |s{ a 



(32) 



(33) 



so that 



■ off 



2^(0) 



?123 



(34) 



In this expression, the meaning of £ 7 (0) is the following. As we said above, the local bias model in Eq. ( |27| ) 
holds for smoothed fields, so ^ (0) is the rms density fluctuation at the smoothing scale. The smoothing filter is 
the Fourier transform, so in this case we should replace ^(0) « A 1 (k) for the average scale k sa (ki+k,2+k3)/3 
under consideration. Using Eq. ( |l8| ) and Eq. (^0|) we find 



(3 



V3tt 



) b 2 



1.53 b 2 



(35) 



Thus, for x 2 initial conditions, the usual Gaussian initial conditions biasing formula is recovered (with no 
additional scale or configuration dependence) provided a proper redefinition of the non-linear bias parameter 
is made. Note that for other spectral indices than n = — 1, the resulting b 2 remains independent of scale. 
In principle, there could be some residual dependence on triangle configuration; however, for n = —1.8 the 
2LPT results described above give negligible residual configuration dependence as well. 



5. Redshift Distortions 

In redshift space, the radial coordinate s of a galaxy is given by its observed radial velocity, a combination 
of its Hubble flow plus "distortions" due to peculiar velocities. The mapping from real-space position x to 
redshift space is given by: 



s = x- f u z (x)z, 



(36) 



where / = d In D\ jd In a w 57 6 , and u(x) = —v(x)/(Hf), where v(x) is the peculiar velocity field, and 
H(t) = (l/a)(da/dT) = Ha is the conformal Hubble parameter (with FRW scale factor a(r) and conformal 
time r). In Eq. (Bq), we have assumed the "plane-parallel" approximation, so that the line-of-sight is taken 
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as a fixed direction, denoted by z. Using this mapping, the Fourier transform of the density field contrast in 
redshift space reads (Scoccimarro et al. 1999) 



S s (k) 



d 3 x 



_ e -ik-x e ifkzU z (x) 



S(x) + f\7 z u z (x) 



(37) 



This equation describes the fully non-linear density field in redshift space in the plane-parallel approximation. 
In linear perturbation theory, the exponential factor becomes unity, and we recover the well known formula 
(Kaiser 1987) 



S s (k) = S(k) (1 + /M 2 ), 



(38) 



where fi = k z /k. Redshift distortions are trivial to include for Q{ 2 3; since only linear PT is involved. 
Assuming linear biasing, Eq. (p8|) gives the monopole of the reduced bispectrum 



n i (l + /?Mi)(l+/3/i|)(l + /?Ml) v Qi2 3 

Vs 123 — ~ ~ <~ X 



(1 + 2/3/3 + /3 2 /5) 2 



A* 



?123 



(39) 



where [3 = f/b±, /Ltjfcj = fej • z and the bar denotes angular average over triangle orientations. Note that 
this results holds irrespective of the type of non-Gaussian initial conditions. After some algebra, we obtain 
A s = C a /(1 + 2/3/3 + /? 2 /5) 2 with 



C, 



/3 2 (7 + 3/3) 
210 



i-6 _i_ 1,6 _i_ jL6 _ L.4L.2 



h 2 h 4 ' — h 4 h 2 

1 2 1 3 



h A h 2 
"-2 ""3 



h 2 h 4 



ft, 2 ft-3 



h 2 h 2 h 2 



(40) 



In Fig. ||we show the correction factor for redshift distortions v4 s . Unlike the Gaussian initial conditions 
case, where the redshift-space correction is only about 10% (Hivon et al. 1995), for non-Gaussian models it 
can be significantly more, depending on the value of /3. The reason is simply due to the scaling of the pri- 
mordial bispectrum being dimensional rather than hierarchical (as that generated by gravity from Gaussian 
initial conditions). That is, from Eq. (|39| ) one expects that approximately A s w 1/^/1 + 2/3/3 + /3 2 /5, which 
describes very well the mean value of A s shown in Fig. [| 

As seen from Eq. (^), the density field in redshift-space is exponentially sensitive to the velocity field. 
Thus, expanding this effect by linear PT is only valid at very large scales. In order to incorporate the 
non-linear effects of the redshift-space mapping, we ran 2LPT realizations in redshift space, which take 
into account this mapping exactly. We assume a cosmological model with £1 = 0.3 and Cl\ = 0.7, for 
which (3 m 0.51. Figures ^ and |5| show the results corresponding to the same triangles as in Figs. | and |, 
respectively. The amplitudes in redshift space are changed, but the overall behavior is the same: the \ 2 
model shows larger bispectrum amplitude and scale dependence than the Gaussian initial conditions case. 
Note that 2LPT calculations in the latter case also automatically include non-linear effects besides those 
due to redshift-space distortions ("loop corrections"), which were not present in the predictions shown as 
dot-dashed lines in Figs. [I] and ^ this explains the mild scale dependence of the equilateral bispectrum in 
Fig. ^| opossite to that of the \ 2 model. In summary, we conclude that the results of Section 3 are robust 
against the effects of redshift distortions. 



Another signature of primordial non-Gaussianity is provided by the rcdshift-space power spectrum. As 
we discussed above, when initial conditions are non-Gaussian, new couplings become available and non-linear 
corrections become important at larger scales than in the Gaussian initial conditions case. In redshift-space, 
this manifests itself in a particular way, as we now discuss. Consider the redshift-space power spectrum, 
from Eq. (|37]) we can write a simple expression for the power spectrum in redshift space (Scoccimarro et al. 
1999) 



d 3 r 
(2tt) 3 



8{x) + fV z u z {x)] [S(x') + /V;«,(sc')]), (41) 



where r = x — x' . This is a fully non-linear expression, no approximation has been made except for the 
plane-parallel approximation. The factors in square brackets denote the amplification of the power spectrum 
in redshift space due to infall, and they constitute the only contribution in linear PT, giving Kaiser's (1987) 
result 

P s (k, fJ ,) = (l + f» 2 ) 2 P(k). (42) 

The anisotropy of the power spectrum in redshift space can be characterized by expanding it in multipole 
moments with respect to /x, the cosine of the angle between a given wave vector fe and the line of sight z. 
This gives a monopole P (k) = (1 + 2/3/ + / 2 /5)P(fc), and quadrupole P 2 (k) = (4/3/ + 4/7/ 2 )P(/c). The 
quadrupole to monopole ratio P2/P0 is thus sensitive to the matter density parameter £1, in fact, assuming 
linear bias Pz/Pq is just a function of p = f/b% (Kaiser 1987; Hamilton 1992). 

On the other hand, at smaller scales the pairwise velocity along the line of sight in the exponential factor 
in Eq. ( |4l| ) starts to play a role. This eventually leads to a decrease in monopole and quadrupole power 
with respect to the linear contribution; in particular, the quadrupole changes sign and becomes negative. 
Note that this effect is distinct from that of velocity dispersion associated with clusters which takes place 
at smaller scales than considered here, and leads to a strong negative quadrupole. In fact, at the scales we 
work the decrease in the quadrupole to monopole ratio can be understood from perturbative dynamics, as 
has been noted before (Taylor & Hamilton 1996; Fisher & Nusser 1996; Hui, Kofman & Shandarin 1999; 
Scoccimarro et al. 1999). 



From Eq. (41), one expects that the scale of quadrupole zero crossing is very sensitive to the magnitude 
of infall velocities. In redshift-space, infall into large-scale structures can become large enough that "near" 
and "far" sides of large-scale structures in real space reverse sides in redshift-space (Hui et al. 1999). At this 
point statistical isotropy is recovered and thus the quadrupole vanishes. The magnitude of infall at a given 
scale k is essentially given by the power spectrum Pg(fc) of the velocity divergence, (9(x) = V • it), which 
evolves according to (compare with Eq.^) 



P e (k) = P J (fc) + 2 j d?qG 2 {k + q,-q) B\k,q), (43) 

where the kernel G 2 (ki, k 2 ) = 3/7+l/2cosd(ki/k 2 + k 2 /ki) + 4:/7 cos 2 9 describes the second-order evolution 
of the velocity field, similar to F 2 in Eq. (]l|) for the density field. Thus, the velocity divergence power spectrum 
is sensitive to the primordial bispectrum. For the \ 2 model with n = — 1 model as in Section 3 gives (compare 
with Eq|J|) 
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Thus, for x 2 initial conditions we expect large-scale infall velocities to be smaller than for Gaussian initial 
conditions. The effect of primordial non-Gaussianity on the pairwise velocity distribution has been explored 
by Catelan & Scherrer (1995). 

Figures || and |?] show results of 100 2LPT realizations of Gaussian and x 2 initial conditions with spec- 
tral index n = —1.4, k n i = 0.33 h/Mpc, £1 — 0.3 and J7a = 0.7. The dashed lines show the predictions of 
linear PT, Eq. (fl2|). Clearly, non-linear effects are important even at scales as large as k = 0.06 h/Mpc. 
Comparing the Gaussian (squares) and the \ 2 model (triangles), we see that the x 2 model shows a smaller 
amplification and the zero crossing of the quadrupole happens at smaller scales, consistent with infall ve- 
locities being smaller as described above. In fact, the quadrupole has not yet reached its large-scale limit 
at k = 0.06h/Mpc. We have checked that this result is not an artifact by verifying that lowering the power 
spectrum normalization brings the quadrupole to agreement with linear PT. Therefore, the shape of the 
quadrupole to monopole ratio is sensitive to primordial non-Gaussianity in this case becoming flatter than 
the corresponding Gaussian model with the same power spectrum. At small scales sensitive to virial motions 
inside clusters, however, we expect a different behavior. Models with primordial positive skewness such as 
the x 2 model tend to have more prominent "fingers of god" (Weinberg & Cole 1992). 

6. Conclusions 

We have explored the predictions of x 2 initial conditions for clustering statistics, in particular the power 
spectrum and bispectrum. We found that extrapolation of the initial conditions using linear perturbation 
theory is not accurate enough at the scales best probed by current and future galaxy redshift surveys 
(k > 0.05 h/Mpc). This is not surprising, non-Gaussian initial conditions provide additional non-linear 
couplings otherwise forbidden by symmetry, thus non-linear gravitational corrections can become important 
at larger scales than in the Gaussian case. 

We showed that when non-linear corrections are included, the bispectrum shows a similar configuration 
dependence as in the Gaussian case; however its amplitude is much larger than the latter and has a residual 
scale dependence not present in evolution from Gaussian initial conditions. We included the effects of galaxy 
biasing and showed that the galaxy bispectrum obeys the same relation to the dark matter bispectrum as 
in the Gaussian case with a proper redefinition of the non-linear bias parameter. Thus, a large linear bias 
decreases the dependence of the bispectrum on triangle configuration, and a non-linear bias contributes a 
constant independent of configuration just as in the Gaussian case. The effects of redshift distortions were 
shown to change the overall amplitude of Gaussian and non-Gaussian bispectra, but clear difference between 
them remains. Thus, we conclude that the bispectrum is a very useful statistic to probe the Gaussianity of 
initial conditions, at least for models where the scaling is dimensional as in the x 2 case. Application of these 
results to current galaxy surveys will be reported elsewhere (Scoccimarro et al. 2000). 

We also discussed the effects of non-Gaussianity on the power spectrum in redshift space. The shape 
of the quadrupole to monopole ratio is very sensitive to the large-scale pairwise velocity which depends on 
the primordial bispectrum through non-linear corrections. For the x 2 model, this leads to a smaller infall 
velocity and thus a supression of overall power and smaller quadrupole zero crossing scale compared to the 
case of Gaussian initial conditions. This provides an additional signature of primordial non-Gaussianity that 
can be explored with the future generation of galaxy redshift surveys. 

I thank Francis Bernardeau, Josh Frieman, Enrique Gaztanaga, Lam Hui, Lev Kofman, and Jim Peebles 
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Fig. 1. — The reduced bispectrum Q for triangles with sides ki = 0.068 h/Mpc and k^ = 2ki as a function 
of the angle 9 between fci and &2. Triangles denote linear extrapolation from x 2 initial conditions, whereas 
square symbols show the result of non-linear evolution. Dot-dashed lines show the predictions of non-linear 
PT from Gaussian initial conditions with the same initial power spectrum as the x 2 model. 
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Fig. 2. — The reduced bispectrum for equilateral triangles as a function of scale k. Line styles as in Fig. |l| 
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Fig. 3. — The correction factor for redshift distortions to the primordial bispectrum, A s , as a function of 
angle 6 between fci and k 2 for triangles with sides fci/fe = 2, for different values of /3 w f2°- 6 /£>i. 
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Fig. 4. — Same as Fig. [I] but in redshift space. 
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Fig. 5. — Same as Fig. ||but in redshift space. 
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Fig. 6. — The power spectrum monopole (top set of curves) and quadrupole (bottom) as a function of scale 
for scale-free initial spectra n = — 1.4 for Gaussian (squares) and x 2 initial conditions (triangles). The dashed 
lines show the predictions of linear PT. 
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Fig. 7. — The power spectrum quadrupole to monopole ratio A2/A0 as a function of scale for scale-free 
initial spectra n = —1.4 for Gaussian (squares) and \ 2 initial conditions (triangles). The horizonal line 
denotes the large-scale limit expected in linear perturbation theory. 



